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Abstract 

Persistence of stationary and traveling single-humped locahzed solutions in the spatial discretiza- 
tions of the nonlinear Schrodinger (NLS) equation is addressed. The discrete NLS equation with the 
most general cubic polynomial function is considered. Constraints on the nonlinear function are 
found from the condition that the second-order difference equation for stationary solutions can be 
reduced to the first-order difference map. The discrete NLS equation with such an exceptional non- 
linear function is shown to have a conserved momentum but admits no standard Hamiltonian struc- 
ture. It is proved that the reduction to the first-order difference map gives a sufficient condition for 
existence of translationally invariant single-humped stationary solutions and a necessary condition 
for existence of single-humped traveling solutions. Other constraints on the nonlinear function are 
found from the condition that the differential advance-delay equation for traveling solutions admits 
a reduction to an integrable normal form given by a third-order differential equation. This reduction 
also gives a necessary condition for existence of single-humped traveling solutions. The nonlinear 
function which admits both reductions defines a two-parameter family of discrete NLS equations 
which generalizes the integrable Ablowitz-Ladik lattice. 



1 Introduction 

We address spatial discretizations of the nonlinear Schrodinger (NLS) equation in one dimension, 

^■"f + + 2|upn = 0, x e M, teM (1.1) 
which has a family of traveling wave solutions 

u = ^ sech(^(x - let - s)) e^^(^-^*)+^^*+^«, (L2) 

where lo G ]R_|_ and (c, s, 6*) € are free parameters of the solution family. The discrete counterpart of 
the NLS equation takes the form: 



h? 

where /i^ > is parameter, and / : C'^ C is a non-analytic function with the properties: 



m„ + -^^±^^ ^ — + = 0, n E Z, t e M, (1.3) 
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PI (symmetry) f{v,u,w) = f{w,u,v) 
P2 (continuity) f{u,u,u) = 2\u\'^u 

P3 (gauge covariance) /(e*"f , e*°ti, e^'^w) = e^'^f{v, u, w) Vq G M 
P4 (reversibility) f{v,u,w) = f{v,u,iv) 

We refer to the model (ll.3t as the discrete NLS equation or simply f/ie A^LS lattice. The discrete NLS 
equation (ll.3t is a symplectic semi-discretization of the continuous NLS equation (ll.lt . where the sec- 
ond partial derivative is replaced with the second-order central difference on the grid x = nh, n £ Z and 
the nonlinearity incorporates the effects of on-site and adjacent-site couplings IOWW04I . From another 
point of view, the discrete model (II. 3t is derived in various branches of physics, most recently, in the 
context of Bose-Einstein condensates in periodic optical lattices LKREOlil . 

The symmetry property PI is required to ensure that the semi-discretization is symplectic. The conti- 
nuity property P2 is needed if the discrete NLS equation (II. 3t is to be matched to the continuous NLS 
equation dl.ll i in the limit h 0. The gauge covariance and reversibility properties P3-P4 originate 
from applications of the discrete and continuous NLS equations to the modelling of the envelope of 
modulated nonlinear dispersive waves in a non-dissipative system IKRBOll . We enforce two additional 
properties on the nonlinear function /: 

P5 f{v, u, w) is independent on h 

P6 f{v, u, w) is a homogeneous cubic polynomial in {v, u, w) 

Neither analysis nor applications require properties P5-P6. However, we use them to limit the search 
of all possible exceptional NLS lattices to a finite-dimensional parameter space. Indeed, we find imme- 
diately that all properties P1-P6 are satisfied if and only if the nonhnear function /(n„_i, u^, Un-i) is 
represented by the ten-parameter family of cubic polynomials, 

/ = ai\Un\'^Un + 0:2\Un\'^{Un+l + Un-l) + a3ul{Un+l + Un-l) + a4(|^in+lP + \Un-l\'^)Un 



+as{\Un+l\'^Un+l + |n„_i|^n„_i) + ag{ui_^^iUn-l + Un+iul-l) 

+aio(|n„+ipn„__i + |-u„„ipn„+i), 



'n 



(1.4) 



where the real- valued parameters (ai, aio) satisfy the continuity constraint: 



«! + 2q2 + 2a3 + 2a4 + 2a5 + 2aQ + ay + 2as + 209 + 2aio 



= 2. 



(1.5) 



This family generalizes the cubic on-site lattice when 



(dNLS) f = 2\un\^u.. 



'n 



(1.6) 
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and the integrable Ablowitz-Ladik (AL) lattice 



(AL) / = \Un\'^{Un+l + Un-l). (1-7) 



A more general nonlinear function of the form (ll.4t was derived in the context of modeling of the 
Fermi-Pasta-Ulam (FPU) problem (see Eqs. (1 1)-(12) in ICKKS93I ) with 

1 1 

oi = 02 = a4 = —, as = ttQ = as = — , a^ = aj = ag = aio = 0. (l-o) 
4 8 

We shall consider persistence of traveling wave solutions (ll.2t in the discrete NLS equation (ll.3t for 
sufficiently small h depending on the nonlinear function (ll.4t . Since both the translational and Gallileo 
invariances are broken, existence of a stationary solution for some u; G Ii C M, 

Unit) = (t){hn)e''^\ (/> : Z ^ C, (1.9) 

does not guarantee existence of a traveling solution for some (w, c) G /2 C M^, 

Unit) = (pihn - 2ct)e''^\ (p:R^C, (1.10) 



where lo is frequency and c is velocity of traveling solutions. In both cases il.9i and (ll.lOt . we are 
only interested in existence of single-humped localized solutions 0(2), which correspond to the sech- 
solutions (II. 2t in the limit h ^ 0. In what follows, stationary and traveling solutions stand for single- 
humped localized solutions unless the opposite is explicitly specified. Direct substitutions of (ll.9t and 
(II. lot to the discrete NLS equation (ll.3t show that the sequence {0n}nez with (j)n = (pihn) satisfies 
the second-order difference equation, 

w</)„ + /(0„_i,0„,0„+i) = 0, n E Z, (1.11) 

while the function (piz) with z = hn — 2ct satisfies the differential advance-delay equation, 

2ic^'iz) = '^^^ + ~ ^{['^ + - o^^iz) + fi^iz - K),4>iz),4>iz + h)), zGM. (1.12) 



Existence of stationary solutions of the second-order difference equation dl.llt can be established for a 
large class of nonlinearities (ll.4t by using the variational method IP05I . Two single-humped solutions 
exist in a general case: one solution is symmetric about a selected lattice node (say n = 0) and the other 
solution is symmetric about a midpoint between two adjacent nodes. 

There are two obstacles for the stationary solution ( ll.9t to persist as the traveling solution ( ll.lOt . First, 



the solution of the stationary problem dl.llt corresponds generally to a piecewise continuous solution 
(j)iz) of the traveling problem (I1.12t on z G M IAEHV05L Jump discontinuities of (piz) may occur in 
between two adjacent nodes, e.g. between Zn = hn and Zn+i = (n + l)h. Second, even if there exists 



a continuous solution ^(z) of the traveling problem (I1.12t with c = 0, this solution may not persist as 
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a continuously dijferentiable solution of the problem (I1.12t with c 7^ 0. The stationary solution 
of the second-order difference equation (ll.llt is said to be translationally invariant if the function 
4>n = 4'{nh) on n G Z can be extended to a one-parameter family of continuous solutions (f){z — s) 
on 2: G M of the advance-delay equation (I1.12t with c = 0, where s G M is an arbitrary translation 
parameter. 

The persistence problem for traveling waves in discrete lattices was recently considered in a number 
of publications. Kevrekidis |K03| suggested a method to define an exceptional nonlinear function 
f{un-i,UmUn+i) from the condition that the discrete NLS equation (II. 3t preserves the momentum 
invariant, 



No characterization of the nonlinear functions (ll.4t that preserve the momentum conservation (I1.13t 
was made in IK03I except for the integrable AL lattice (ll.7t . In addition, it was not shown in IIK03I 
that the preservation of momentum (I1.13t guarantees the existence of a continuous function (/){z) in the 
solution (II. lot with c = 0. 

Oster, Johansson, and Eriksson IQJE03I addressed the most general Hamiltonian discrete NLS equation 
in the form. 



where F : M is a symmetric potential function of quartic polynomials which satisfies the gauge 

invariance property: F(e*"u, e^'^w) = F{u, w) Ma G M. The dNLS lattice (fTel and the FPU lattice (fTSl 
are particular examples of the Hamiltonian discrete NLS equation ( I1.14t . It was found in IIOJE0 3 1 that 
there exists a configuration in parameters of quartic polynomials F{un, Un+i) which provides a reduced 
Peierls-Nabarro potential and enhanced mobility of localized modes. No analysis of the locahzed modes 
with enhanced mobility was developed in IOJE03 I. 

Ablowitz and Musslimani IAM03II outlines problems in the numerical approximations for traveling 
wave solutions dl.lOl) with c 7^ in the context of earlier contradictory publications (see also review in 
IPR05I ). An asymptotic method was developed in IAM03I to show that if the solution exists for c = 
it can be continued to non-zero values of c as perturbation series expansions in powers of c. 

Pelinovsky and Rothos IIPR05I computed the normal form for bifurcations of traveling solutions dl.lOt 
near the special value of the parameter c = 1/h. The normal form is represented by the third-order 
ODE related to the third-order derivative NLS equation IIPR05I . Existence of embedded solitons in 
the third-order derivative NLS equation was considered in the past (see IIYA03I IPY05I and reference 
therein) and the previous results suggested that the dNLS lattice ( ll.6t did not support bifurcation of 
single-humped traveling solutions in a sharp contrast to the case of the AL lattice (ll.7t . {Double-humped 
traveling solutions in the normal form reduction were predicted for the dNLS lattice M.b\ in IIPR05I and 




(1.13) 




dH 




(1.14) 
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numerical approximations of the double-humped traveling solutions in the differential advance-delay 
equation (I1.12t were obtained in IC06I .') 

We shall review the results and conjectures of the previous works IAM03I IK03I IQJE03I IPR05I in the 
context of the cubic nonlinear function ilAi . Our results are based on the recent works |BOP051 and 
IDKY05I . where similar ideas are developed for monotonic kinks of the discrete c/)'^ theory. See also 
IOPBO5IIIPO6II for other relevant results on monotonic kinks in discrete Klein-Gordon lattices. 

In Section 2, the class of exceptional nonlinear functions in the general family of cubic polynomials (ll.4t 
is identified from the condition that the second-order difference equation dl.llt admits a reduction to 
the first-order difference equation. It is shown that the discrete NLS equation (II. 3t with an exceptional 
nonlinearity preserves the momentum conservation (I1.13t but admits no Hamiltonian in the form (11.141) . 

In Section 3, it is proved that the first-order difference equation admits a one-parameter family of station- 
ary solutions ( 11.91 ) for any u; > and sufficiently small h, and these stationary solutions are translation- 
ally invariant. Unlike the case of monotonic kinks in IBOP05I. the proof of existence of a translationally 
invariant single-humped localized solution is based on a construction of two sequences which depend 
continuously on the initial value. One sequence is monotonically decreasing to zero and the other se- 
quence is monotonically increasing until the turning point and it decreases monotonically beyond the 
turning point. The main outcome of this analysis is a clear evidence that the reduction to the first- 
order difference equation gives a sufficient condition for existence of translationally invariant stationary 
solutions and a necessary condition for existence of traveling solutions near c = for tj > 0. 

Section 4 tests the exceptional nonlinear functions by the reduction of the differential advance-delay 
equation ( 11.121 to the normal form near the special values of parameters cj = (vr — 2) //i^ and c = 
1/h. Reductions of the normal form to integrable Hirota IH73II and Sasa-Satsuma IISS91I equations 
and existence of single-humped solutions in the integrable normal form give a necessary condition for 
existence of traveling solutions in the differential advance-delay equation (I1.12t near the special values 
of parameters (w, c). 

Section 5 presents the explicit form of the two-parameter exceptional nonhnear function (ll.4t that passes 
through both tests and may contain families of traveling solutions. This function generalizes the AL 
lattice il.li which is known to admit exact traveling solutions between the two limits c = and c = l/h. 
Numerical test of persistence of traveling solutions in the differential advance-delay equation ( 11.121 ) is 
outhned as an open problem. 

2 Reductions to the first-order difference equation 

We shall consider the stationary solution il.9l of the discrete NLS equation ( I1.3I ). which satisfies the 
second-order difference equation dl.l II ). We obtain the constraints on the function /(0n-i> 4'n, 0n+i) 
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in the family of cubic polynomials ilAl from the condition that the second-order difference equation 
(ll.llt admits a conserved quantity 

E = ■^\<pn+l -<t>n\'^ - {4'n4>n+l + <^n0n+l) + g{4>n, <An+l), (2-1) 

where (7 : i— > M is a non-analytic function and E M is constant in n e Z. Due to properties P1-P6 
for the function /(</)„_i, (pn, 4>n+i), the nonhnear function g{(pn, fpn+i) satisfies the properties: 

51 (symmetry) g{u,w) = g{w,u) 

52 (continuity) g{u,u) = and u) = ^(n, n) = |npu 

53 (gauge variance) g{e^°'u, e^°'w) = g{u, tf) Va G M 

54 (reversibility) g{u, w) = g{u, w) 

55 g{u, w) is independent on h 

56 g{u, w) is a homogeneous quartic polynomial in (n, w) 

The most general polynomial g{4'n, (l>n+i) that satisfies properties S1-S6 takes the form: 

g = ll{\(pn\'^ + \<Pn+l\'^){(l>n+l(f>n + 4>n+l4>n) +72\4>n\'^\4>n+l\'^ 

+73i&n+l + 4>l4>l+l)+7i{\(Pn\^ + I'/'n+ll'), (2.2) 

where the real- valued parameters (71, 72, 73, 74) satisfy the continuity constraint: 

471 + 72 + 273 + 274 = 1. (2.3) 

In the limit /i — > 0, the second-order difference equation (ll.llt reduces to the second-order ODE: 

^" -uj(j) + 2\(p\'^(l) = 0, (2.4) 

while the first-order difference equation (I2.lt transforms to the first integral of (I2.4t 

E=\^'\^-uj\(l)\^ + \^\\ (2.5) 

The first-order ODE (12. 5t with E = admits a unique locahzed solution 4>{z) = ^JZ) sech(y^(2: — s)) 
with s G M and w > 0, which agrees with the exact solution (ll.2t for 2; = x and c = 0. Therefore, 
the first-order difference equation ( 12.11 ) is a symplectic discretization of the first-order ODE ( 12.51 ). This 
precise methodology was used in jpKY05 1 to construct translationally invariant monotonic kinks in the 
discrete (j)^ equation. (A more general approach is reported in IBOP05I .) We first specify on the corre- 
spondence between the difference equations dl.l It and ( I2.lt and then prove three elementary results on 
existence of the reduction (I2.lt with the nonlinear function (I2.2t and its relation to conservation of M 
in (fTBl and H in (fTTO. 
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Lemma 2.1 Let the function g{(pn, fpn+i) with properties S1-S6 be related to the function f {(pn- 1 , <?^'n , (pn+i ] 
with properties P1—P6 by 

1 r 



9i<f)n, 4>n+l) - g{(pn-l,(pn) = ^ {(Pn+l " '/'n-l )/(0n-l , 0n , (j^n+l) + i4>n+l " 0n-l)/(0n-l , (/"n , 4>n+l] 

If the sequence {(pn}nez i^ ciny solution of the second-order difference equation il.llt . it satisfies the 
first-order difference equation ( I2.il) . If the sequence {(j)n}n& i^ ci non-constant real-valued solution of 
the first-order difference equation ( 12. it , it also satisfies the second-order difference equation M.ll\ . 

Proof. By subtracting the first-order difference equation (I2.lt with n = n from that with n = n — 1, 
we have 

-((pn+l - (Pn-l) I 

The first statement follows immediately from the relation between (/>„, (pn+i) and fif((/)„, 0n+i)- 

Let us now assume that the first-order equation (I2.lt is satisfied and functions /{(pn-i, (pn, (f'n+i) and 
9{4'n, (f'n+i) are related as above. Then, the non-constant sequence {(j)n}nez satisfies the second-order 
equation, 

^^2" ^"^'^ ~ + /(0n-l, 4>n-, (pn+l) = i{4>n+l " (t>n-l)H4'n-l, 4>n, 4>n+l), n G Z, 

where /i : C"^ i-^ M is arbitrary function. If the non-constant sequence is also real-valued, then h = 0. U 

Remark 2.2 The one-to-one correspondence between the second-order and first-order difference equa- 
tions dl.l It and ( 12.11 ) exists only for non-constant real-valued solutions. Complex-valued solutions of 
the first-order equation (I2.lt can give solutions of the second-order equation with an additional function 
h{(j)n-i, 4>n, 4>n+i)- This property can be used for a full time-space discretization of the NLS equation 
dl.lt . when the derivative term (l)'{z) in the differential advance-delay equation ( 11.121 ) is replaced by the 
difference term — 0n-i)/(2/i) such that the function h{(j)n-i, ipn, 4'n+i) = c/his constant. We 

will not consider full time-space discretizations of the NLS equation dl.lt in this paper. 

Lemma 2.3 The reduction of the second-order difference equation M.ll\ to the first-order difference 
equation (I2.il ) exists provided that the nonlinear function M.4\ satisfies the constraints: 

= ai — Qfg, Q!5 = CKQ, a-j = OL\ — Iol^^ aio = — ag, (2.7) 

such that (ai, eta, ae, asi ctg) £ are free parameters. 

Proof. By using symbolic computations with Wolfram's Mathematica we find from (Ol, (IZli and 
d2.6t that a\ = 271, 02 = 72> as = 273, and ag = 74 under the constraints d2.7l ). ■ 
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Corollary 2.4 When (f>n € M, the first-order difference equation ( 12.71) is characterized by the symmetric 
quartic polynomial function 

g = Pi(j)l(t)l^^ + P2{4'1 + 0^+i)0n(/>n+l + P3{4>t, + </'n+l)> (2-8) 

where /3i = 72 + 273 = 02 + 03, 02 = 271 = ai, and /^s = 74 = as under the continuity constraint 

(ii + 2/?2 + 2/33 = 1. (2.9) 



Remark 2.5 Corollary l2.4l recovers the result of IBOP05I . where real- valued nonlinear functions of the 
difference equation (II. lU are considered in the context of kink solutions of the discrete 0^ model. 

Lemma 2.6 The discrete NLS equation l\1.3i conserves the momentum invariant ( 17.731) provided that 
the constraints \2.7l are met. 

Proof. Computing time-derivative of (11.131) and using the stroboscopical summation in n, we convert 
the result to the following irreducible remainder: 

M = {ai - - aQ)^^{\Un+l\^ - \Un\^){Un+lUn + Un+lUn)) 

+ (04 - 05 - a^) ^^{\Un+l\^{UnUn-l + UnUn-l) - lUn-l]"^ {UnUn+1 + UnUn+l)) 

+ (as - ae) 

+ ("8 - ag - aio) ^(kn+lP - \Un-l\'^){Un+lUn~l + Un+lUn~l)- 

The constraints ( 12.71 ) give the most general solution of the system of homogeneous linear equations, 
which follows from the momentum conservation M = 0. ■ 

Corollary 2.7 There exists a one-to-one correspondence between the set of nonlinear functions ( 17.41 1 
that conserves the momentum ( 17.731 1 and the set of nonlinear functions that supports the reduction to the 
first-order difference equation ( 12.71 1. 

Remark 2.8 Corollarv l2.7l proves the conjecture of IK03I for the case of a discrete NLS equation (ll.3t 
with the cubic polynomial function ilAl . A similar result for a discrete Klein-Gordon equation was 
obtained in the most general case in IK031 Appendix] by explicit computations. 

Lemma 2.9 The discrete NLS equation ( 17.31 1 has the Hamiltonian structure M.14\ provided that 

a2 = 2a3 = 2a%, 05 = 07 = ag = aio = 0, (2.10) 
such that (ai, 02, 04, ae) ^ are free parameters. 
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Proof. The most general symmetric real-valued gauge-invariant function F{un,Un+i) of the quartic 
polynomials in the Hamiltonian (I1.14t takes the form: 

F = (5i(|n„|^ + |u„+i|^)(u„u„+i + u„n„+i) + 52|w„p|-u„+ip 

Due to the symplectic structure (I1.14I) . the potential function F{un,Un+i) produces the nonlinear func- 
tion f{un-i,Un,Un+i) in the form, 

f = Si [\un+i\'^Un+i + + 2|ii„|^(n„+i + ii„_i) + n^(nri+i + Un-i)] 

Therefore, 02 = 26i, = 62, cxq = 283 and ai = 4J4 under the constraints (I2.10t . ■ 

Corollary 2.10 None of the discrete NLS equations with the Hamiltonian structure M.14\ is momentum- 
preserving, including the dNLS lattice M.61 and the FPU lattice M.8\ . 

Remark 2.11 Corollarv ITTOl agrees with the conclusion of IDKY05I obtained for the discrete c^^ 
model. 



Remark 2.12 The three-parameter discrete NLS equation with the Hamiltonian structure (I1.14t was 
derived in IOJE03I for modeling of arrays of optical waveguides. The model of I OJ E03II corresponds 
to the reduction of the symmetric potential function u„_|_i) with 82 = 4^3. Although it is claimed 

in IIOJE03I from results of numerical simulations that traveling solutions may have enhanced mobility 
in this model, no momentum-preservation is supported and the second-order difference equation (II. lit 
provides no reduction to the first-order difference equation (I2.lt . 

3 Existence of translationally invariant stationary solutions 

We shall consider existence of continuous solutions of the advance-delay equation (11.121) with c = 
from existence of the translationally invariant solutions of the second-order difference equation dl.lll) . 
We show that there exist translationally invariant stationary solutions if the second-order difference 
equation dl.llt is reduced to the first-order difference equation (I2.lt . Thus, the existence of the reduction 
to the first-order difference equation ( 12.11 ) gives the sufficient condition for existence of translationally 
invariant stationary solutions. For simplicity, the analysis is developed for real-valued solutions of 
the first-order difference equation (I2.lt . By Lemma 12.11 such non-constant solutions correspond to 
real-valued solutions of the second-order difference equation dl.lll ). The initial-value problem for the 
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first-order difference equation il.li with = in the space of real-valued localized solutions takes the 
implicit form, 

{(pn+l - (t>nY = h'^L0(j)n4>n+l " h'^ g{(pn, (pn+l) , U £ Z, 

00 = ^, 

where 9? G M is the initial data, iterations in both positive and negative directions of n G Z are con- 
sidered, and 51 : M is given by (I2.8t . Let x = (pn and y = (pn+i and rewrite g{(j)n,4>n+i) 
as 

g{x, y) = pix^y^ + p2xy{x'^ + y^) + /?3(x^ + y^)- 

Due to the continuity constraint (I2.9t . we have g{x, x) = x^ and dxg{x, x) = dyg{x, x) = 2x^. We first 
establish the sufficient condition that the second-order difference equation (ll.llt has only real-valued 
solutions and then develop analysis of the initial-value problem (I3.lt in n G Z. Our main result is 
Proposition B. lOl which guarantees existence of a particular single-humped localized sequence {0n}nez 
to the initial-value problem (I3.lt with to > and sufficiently small h. The sequence corresponds to the 
translationally invariant solution of the second-order equation (ll.llt which converges as /i ^ to the 
solution (j)s{z) = ^/uj sech.{^/ujz) at the points z = hn with n G Z. 

Lemma 3.1 The second-order difference equation HI. Hi with sufficiently small h admits a conserved 
quantity 

J = 0n0n+l - 0n0n+l (3.2) 

for bounded sequences {(j)n}nez if 

ay = 2aQ, aio = as + ag- (3.3) 



Proof. The nonlinear function (11.41) satisfies the relation 

^n/(0ri-l, 0ri, 0ri+l) ^ 0ri / (0ri- 1 , 0n , 0n+l ) = (^n'/'n+l ^ 0ri0n+l + <t>n4>n-l - 0n0n-l) Qn + R-n 

where 

Qn = {0L2 - a^)\4>n\^ + asd^n-lP + \4>n+l\^) + a9(0n-10n+l + 0n-10n+l) 

+a6(0n0n+l + 0n0n+l + ri0n— l) 

and 

Rn = ("7 - 2a6)(0n0n-10n+l " 0n0n-10n+l) 
+ (aiO - as - ag) [{4>n<Pn~l - 4>n<Pn-l)\4>n+l\'^ + (0n0n+l " <t>n4>n+l)\<t>n-l\'^] ■ 

If ar = 2a6 and aio = ag + ag, then i?„ = 0. If the sequence {4>n}n& is bounded, then 1 + h^Qn > 
and J is constant for all n G Z provided that h is sufficiently small. ■ 
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Corollary 3.2 Under the constraints ( 13.31) . real-valued solutions of the second-order difference equa- 
tion M.lll with J = are the only localized bounded solutions {(j)n}n& with (pQ £ M. 

Remark 3.3 While it is not clear if localized complex- valued solutions of the second-order difference 
equation dl.lll ) may exist when 07 / 2aQ and aio / ag + ag, we observe that the reduction to the first- 
order difference equation (12.11) has six free parameters (see Lemma l23l but the nonlinear function in the 
form ( 12.21 ) has only four parameters. Therefore, two parameters can be used to satisfy the constraints 
i3.3l . Combining the constraints i2Jl and i3.3l . we obtain the most general parameterization of the 
nonlinear function (ll.4t that admits the conserved quantities (12.11) and (13. 2t : 

ai = 4a6, a4 = Sag? = a^, 07 = 2aQ, ag = 0, aio = ag, (3.4) 

where {02, 03, oq, ag) G are free parameters. By Lemma ITTI and Corollary 13.21 all localized non- 
constant solutions of the second-order difference equation dl.l II) under the constraints (13.41) are real- 
valued and the sequence {0n}nez is equivalently found from the initial-value problem ( 13. U . 




Figure 1: Solutions of the quartic equation ( 13.51) with /3i = /?2 = and = 0.5. 

Lemma 3.4 Let Q{x,y) = toxy — g{x,y) and u > is parameter There exists /iq > and L > y/uj 
such that the algebraic (quartic) equation 

(y-xf = h^Qix,y) (3.5) 

defines a convex, simply-connected, closed curve inside the domain < x < L and < y < L for 
< h < liQ. The curve is symmetric about the diagonal y = x and has two intersections with the 
diagonal at y = x = and y = x = ^/uj. 
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Proof. The statement of lemma is illustrated on Figure 1. Symmetry of the curve about the diagonal 
y = X follows from the fact that g{x,y) = g{y,x). Intersections with the diagonal y = x follows 
from the continuity constraint that gives Q{x,x) = x'^{uj — x"^). To prove convexity, we consider the 
behavior of the curve for sufficiently small h in two domains: in the interval < x < and in a small 
neighborhood of the point y = x = ^/uj with the radius tq/i^, where xq = ^/uj — tqIp? and tq > is 
/i-independent. Let y = x + z and consider two branches of the algebraic equation (I3.5t in the implicit 
form 



F±{x,z,h) = hyjQ{x, X + z) = 0. (3.6) 

such that for x > 

F±(x,0,/i) = ^hx^/uj-x'^, a^F±(x,0,/i) = 1 T ■ 

It is clear that F±{x, 0, 0) =0 and ^^^^-(x, 0, 0) = 1, while F±{x, 0, h) is continuously differentiable 
in X and h and dzF±{x, 0, h) is uniformly bounded in x and h on any compact subset of x G [0, xq], 
where xq = ^/uj — r^h?' and ro > is /i-independent. By the Implicit Function Theorem, there 
exists ho > 0, such that the implicit equations ( I3.6t define unique roots z = ±hS±{x, h) in the domain 
< X < Xq and < /i < /iq, where hS±{x, h) are positive, continuously differentiable functions in 
X and h. Therefore, the algebraic equation (I3.5t defines two branches of the curve located above and 
below the diagonal y = x. When /iq is sufficiently small, the two branches are strictly increasing in 
the interval < x < xq- In the limit /iq 0, the two branches converge to the diagonal y = x on 
< X < ^/uj. Derivatives of the algebraic equation (I3.5I) in x are defined for any branch of the curve by 

y [2{y - x) -h^iux - dyg{x,y))\ = 2{y - x) + h^{u}y - dxg{x,y)), (3.7) 
y"[2{y-x)-h^{ojx-dyg{x,y))\ = -2{y' - if + h\2u:y' - ~g) , (3.8) 

where g = d^^g{x, y) — 2y'd'^yg{x, y) — {y'Ydyyg{x, y). It follows from the first derivative ( 13.71 ) that 
y' = —1 at the point y = x = ^/lo for any /i > 0, such that there exists a small neighborhood of the 
point y = X = y/uj where the upper branch of the curve is strictly decreasing. Let i?^ be the upper 
semi-disk centered at y = x = ^/uj with a radius r^h?, where ro > is /i-independent. Then, there 
exists a /i-independent constant C > such that 

2{y-x)-h''[u:x-dyg{x,y))>Ch\ (x,y) G 

It follows from the second derivative ( 13. 8t for sufficiently small h that y" < —Ci/h?'m{x,y) G B'^ 
with Ci > 0. Therefore, the curve is convex in B^. By using the rescaled variables in B^: 

x = y/uj-Xh^, y = y/uj + Yh'^, 

we find that the curvature of the curve in new variables is bounded by Y"{X) < —Ci and therefore, the 
first derivative Y'{X) and so y'{x) may only change by a finite number in (x, y) G Bh. Therefore, there 
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exists < C2 < oo such that y'{x) = C2 at a point x, where x < xq = ^JIo — r^h? . By the first part 
of the proof, the upper branch of the curve is monotonically increasing for < a; < xq- By the second 
part of the proof, it has a single maximum for xq < x < y/uj. Thus, the curve defined by the quartic 
equation ( I3.5t is convex for y > x > (and for < y < x by symmetry). ■ 

Corollary 3.5 Let ipQ be a maximal value of y and ipo be the corresponding value of x on the upper 
branch of the curve defined by the quartic equation \3. 51) . Then, ipQ < \pjj < ipo and 

<^o = V^-§V^/i2 + 0(/i^), i;o = + 0{h^). (3.9) 

o o 

Lemma 3.6 The initial-value problem ( 13.71 1 with uj > and < /i < /iq admits a unique monotonically 
decreasing sequence {(j)n}^=ofor any < ip < ^/uj that converges to zero from above as n ^ 00. The 
sequence {(j)n}'^=Q is continuous with respect to h and (p. 

Proof. By Lemma (see Figure 1), there exists a unique lower branch of the curve in (13. 5t below 
the diagonal y = x for < x < ^/uJ and the monotonically decreasing sequence {4'n}^=o with 
< (f < satisfies the initial-value problem: 

f (t>n+l = (t>n- hS^{<pn, h), 71 G N 

\ (t>0 = f 

where hS- (0, h) > is continuously differentiable with respect to (p and h. We shall prove that the 
monotonically decreasing sequence {</>n}$^o converges to zero from above. Since Q{x,y) is a quartic 
polynomial, there exists a constant C > that depends on lu and is independent of h, such that 



?<Ch^(t>l (3.11) 



for all 4>n+i < (pn- If h is sufficiently small, such that C/i^ < 1, then < 0n+i < 4>n, and the sequence 
{4>n}^=o is bounded from below hy (p = 0. By the Weierstrass Theorem, the monotonically decreasing 
and bounded from below sequence {i?!>„}^q converges as n — > oo to the fixed point = 0. Continuity 
of the sequence {0„}^q in h and (p follows from smoothness of hS- {4>, h) in h and ip. ■ 

Lemma 3.7 There exists N > 1, such that the initial-value problem ( 13. il l with u; > and < h < Iiq 
admits a unique monotonically increasing sequence {(j)n}n=~oof^^ '^^y < < y/uj that converges to 
zero from above as n ^ —oo. The sequence is continuous with respect to h and ip. 

Proof. By Lemma l3!4l (see Figure 1), there exists a unique upper branch of the curve in (13. 5t above 
the diagonal y = x for < x < and the monotonically increasing sequence with 
< < satisfies the initial-value problem: 

f 0„+i = (/'n + /i5'+ h), (-n) G N 
1 (/)0 = 
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where hSj^{4),h) > is continuously differentiable with respect to cf) and h. Existence of N > 1 
follows from the same equation (I3.12I) for < n < — 1. The proof that the monotonically increasing 
sequence {(t)n]n=-oo converges to zero from above as n — > — oo is similar to the proof of Lemma 13.61 
Continuity of the sequence {(l)n}n=-oo ^ ^^'^ "-P follows from smoothness of hSj^{(j)^ h) in h and ip. ■ 

Lemma 3.8 The initial-value problem \3.1\ with w > and < h < ho admits a unique 2-periodic 
orbit {<l)n}n& with cpn+i / 4>n and 0„+2 = 4>nfor any ipo < ip < itQ. 

Proof. By Lemma lT4l (see Figure 1), the curve in (I3.5t is symmetric about y = x and has two branches 
in X below y = x for ^/tJ < x < -00- Therefore, any initial data on the branch between ((/jq, V'o) and 
{tpo, po) leads to a unique 2-periodic orbit. ■ 

Corollary 3.9 The initial-value problem (lj.il) with uj > and < h < ho admits the following 
particular sequences: 

• For any given < (p < ipo, there exists a symmetric single -humped sequence {4>n}n& with 
maximum atn = 0, such that (pn = (^^^ Figure 2(a)). The single-humped sequence is unique 
for ip = ipo (see Figure 2(b)). Let Son denote the countable infinite set of values of {(pn}n&for 
the single-humped solution with p = ipo- 

• For any given ^/uj < p < ipo, there exists a symmetric double-humped sequence {(t)n}ne'L with 
local minimum at n = 1 and maxima at n = and n = 2, such that (pn = <P-n+2 and piQ < 4>i < 
a/w ( see Figure 3(a)). The double-humped sequence becomes a unique 2-site top single-humped 
sequence for p = y/uj (see Figure 3(b)). Let S'off denote the countable infinite set of values of 
{4'n}nezfor the 2-site top single-humped solution with p = ^JZj. 

• For any p G (0, V'o)\{'S'oni ■S'ofr}. there exists a unique non-symmetric single-humped sequence 
{(pn]n& with 4>k / (pmfor all k^m. 

Proposition 3.10 The second-order difference equation U.llt admits a translationally invariant single- 
humped sequence {(t)ri\n& for any < 0o !^ V'o with 0„ = (j){nh — s), where n S Z, s S M, and 
(t){z) is a continuous function. In the limit h ^ 0, the function (j)(^z) converges pointwise to the function 
(psiz) = sech(-^/Li7z), i.e. there exists C > and s G M such that 

max — (psinh — s)\ < Ch^. (3.13) 
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Figure 2: Left: a construction of the symmetric single-humped sequence from the solution of the quartic 
equation (I3.5I) . Right: an example of the symmetric single-humped solution of the difference equation 
(I3.lt for (/9 = -00- 




Figure 3: Left: a construction of the symmetric double-humped sequence from the solution of the quartic 
equation (I3.5t . Right: an example of the 2-site top single-humped solution of the difference equation 
(I3.lt for Lp = y/uj. 
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Proof. By the geometric construction of Corollary 13.91 there exists a single-humped sequence {(j)n}nez 
for any < (po < i^o, which consists of the unique symmetric sequence for (po G ^on, the unique 2-site 
top sequence for (f>o G 5'off and the unique non-symmetric sequence for 0o £ (0, ilJo)\{Son, 5'ofr}- By 
Lemmas l3.6l and im the single-humped sequence is continuous in h and (pQ such that it is translationally 
invariant. By Lemma ITTl the real- valued non-constant sequence satisfies identically the second-order 
equation dl.lll ). It remains to show the pointwise convergence of the sequence {i?i>„}nez to the sequence 
{(psinh — as /i ^ 0, where (j)s{z) = -y/u; sech.{^/uJz). By the translational invariance of 

the sequence, we can place the maximum at n = 0, such that cpo = iJjq, which is equivalent to the 
choice s = 0. Then, the bound ( 13.131 with s = coincides with the bound for the convergence of the 



symmetric single-humped solution with {(pn}nez = 5'on of the second-order difference equation (II. lit 
to the solution (psihn) = ^JTo sech (y^/in) of the second-order ODE MA\ . The standard proof of the 
error bound (I3.13t is performed with a formal power series in which involves hyperbolic functions 
for exponentially decaying solutions (see Appendix B in IBGKM91I ). The formal series is different 



from the rigorous convergent asymptotic solution of the difference map dl.lll ) along the stable and 
unstable manifolds by the exponentially small in h terms (see ITOOI for application of this technique to 
the second-order difference equation with a quadratic nonlinearity). We note that a similar result (but 
with a different technique based on analysis of Fourier transforms) was proved in IFP99I for the discrete 
Fermi-Pasta-Ulam problem. ■ 

Example 3.11 Let us consider the explicit example of the AL lattice (ll.7t . when fi\ = \ and /32 = /^s = 



0. In this case, the second-order difference equation dl.lU admits the exact solution for translationally 
invariant single-humped solutions: 

'Pn = "00 sech(K/in — s), (3.14) 
where (s, k) G are arbitrary and (i/^O) t^) are defined by 

V'o = r — ^ = 79 sinh^ • (3-15) 



h ' h? V 2 , 

The set 5*011 for the symmetric single-humped solution is defined by values of {4>n}n& for s = (see 

nh 
2 



Figure 2(b)) and the set S'off for the 2-site top single-humped solution is defined for s = We find 



that 



^ fKh\ / a;2/j2 
V'o = vwcosh ( — I = + (3.16) 

in agreement with the asymptotic formula d3.9l l. The algebraic equation d3.5t can be reduced to the 
explicit roots for the one-step maps hS± (0, h) in the form: 

= W^hW) ' 
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which satisfy all properties derived in Lemma l34l In particular S^{(/), h) > for < i?i> < ^/uj and 
S'_((/), /i) > has one root for < < y/uj and two roots for ^/uJ < cf) < ^q, where ipQ is given by 
( I3.16I ). In the limit /i ^ 0, the expression ( I3.14I ) with ( I3.15I ) converges to the solution ( I1.2I ) with u) = 
and c = 0. 

Remark 3.12 For real-valued non-constant solutions, the second-order difference equation dl.llt can 
be reduced to another first-order difference equation: 

-((/)n+l - (pn) = 9i4>n, (pn+l), (3.17) 

where 4>n+i) is a symmetric quadratic polynomial. The form ( 13.171 ) was used in IBOP05II in a 
search of another family of exceptional discretizations which admits translationally invariant mono- 
tonic kinks in the discrete (j/^ equation. By the Implicit Function Theorem, the first-order difference 
equation ( 12.11 ) has only one branch of solutions on the plane (f^n+i)- As a result, one can not con- 
struct two (decreasing and increasing) sequences {(/)n}^Q in the first-order difference equation (I3.17t 
for translationally invariant single-humped solutions. Therefore, other exceptional discretizations of 
IBOP05I are irrelevant for localized solutions of the discrete NLS equation ( ll.3t . 



4 Existence of traveling solutions near uj = {tt — 2)/h? and c= 1/h 

We shall consider existence of continuously differentiable solutions of the differential advance-delay 
equation ( 11.121 ) with c ^ 0. A convenient parametrization of the solution is represented by the transfor- 
mation of variables 



and parameters 



(z) = ^$(Z)e^/^^, Z = ^, (4.1) 
n n 



2 2 1 

Lu = — /3c + (cos (3 cosh k — 1) , c = - — sin (3 sinh k, (4.2) 
n /i^ fiK 



such that the pai^ameter h is scaled out the new equation for ^{Z): 
Sinn K 

2i sin (3 + 2 cos /3 cosh k^{Z) 

K 

= ^{Z + l)e*^ + - l)e-'^ + / - l)e-^^, ^>(Z), $(Z + l)e^^) , Z G M. 

If the single-humped locahzed solutions to the differential advance-delay equation exist, then the func- 
tion ^{Z) decays in Z with the real-valued rate k, i.e. ^{Z) ~ e^'^l'^l. (Near the bifurcation line «; = 0, 
the function ^{Z) may also contain decaying exponential terms with oscillatory tails due to presence 
of complex eigenvalues but these terms decay faster than the tail ^{Z) ~ e"*^'^'.) The boundary on 
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possible existence of traveling solutions corresponds to the line k = 0, such that the traveling solutions 
may only exist in the exterior domain to the curve 

io = ^pc- ^sin^ ^, c = -^ sin/3, /3e[0,27r]. (4.3) 

The boundary on existence of traveling solutions is shown on Figure 4 for c > 0. We apply the technique 
of the normal form reduction from IPR05I when parameters {to, c) are close to the special values co = 
{'K — 2)/h? and c = \/h, which correspond to (3 = 7r/2. By using a modified transformation of variables 
and parameters 



we rewrite the differential advance-delay equation (I1.12I) in the equivalent form, 

i ($(C + 6) - $(C - 6) - 26$'(C)) = e- {2W^\C) + ^^{0) - e^i-m - 6), $(C), i$(C + e)), 

where (fi, V) G are rescaled parameters {lu, c) G /2 C M? and properties P3 and P6 of the function 
f{v, u, w) are used. The small parameter e > defines deviation of parameters (w, c) from the values 
w = {'K — 2)/h? and c = as well as the amplitude and localization of the solution 4>{z). By using the 
formal Taylor series expansions of the shift operators <I>((^ it e) in powers of e, the differential advance- 
delay equation for $(C) is converted to a third-order ODE which is related to the third-order derivative 
NLS equation JPROS |. The formal reduction can be proved with the rigorous technique of the center 
manifold and normal forms (see analysis in IIP06I for kinks in the discrete (j)'^ equation). 

The formal reduction of the hnear and nonlinear parts of the differential advance-delay equation for 
<I>(C) leads to the expansions, 

$(C + 6) - <I>(C - 6) - 26<I>'(C) = ^<i>'"(C) + 0{t>), 

and 

f{-i^{Q - e), ^(C), ««'(C + e)) = (ai + 2a4 - 205 - 2a6 + aj) |«>|2$ 
+2ie (02 + 2a8 - 209) - 2ie (03 - ag - ag + aio) ^>^^'(C) + 0{e^). 

By rescaling the amplitude of $(C) one can bring the term |<I>p<I> to the order of e^. However, the 
third-order ODE 

3 ' ' ' 

has no single-humped localized solutions for any (17, V) G (see references in IYA03llPR05l l There- 
fore, the necessary condition for existence of single-humped travehng solutions is the constraint on 
parameters of the cubic polynomial function (ll.4t : 

ai + 2a4 - 2a5 - 2aQ + aj = Q. (4.4) 
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When other constraints (I2.7t and (I3.3I) are taken into account, the new constraint produces 06 = 
in the representation (I3.4t . Under the constraints (I2.7I) . (I3.3I) . and (I4.4I) and the normalization (I1.5I) . 
no rescaling of the amplitude of ^{C) is needed and the truncated normal form for traveling solutions 
becomes 

- 2iV<^>' + 2i|^>|2$' + i7$(|$p)' = 0, (4.5) 

3 

where 7 = —2a^ = 2(02 + 203 — 1) is a real parameter. Existence of single-humped localized solutions 
in the third-order ODE ( I4.5t is related to existence of embedded solitons in the third-order derivative 
NLS equation (see recent survey in IPYOSI l We shall represent the basic facts about existence of single- 
humped localized solutions of the third-order ODE i4.5l . By using the transformation of variables and 
parameters 

$ = A^(Z)e~*^^ z = xc, n = -k(X^ + k^), V = - (X^ - 3k^) , 

3 6 

we rewrite the ODE ( 14. 5t in the form, 

— - ^' + el^I-p^' + 37^(|^'p)') + k - ^' + 2|^|2^) = 0, (4.6) 
3 

where the real-valued parameters (A, k) are arbitrary. The value A = defines the bifurcation cutoff in 
the family of localized solutions of the ODE (14.61) . which is equivalent to the curve V = ythr(^) in the 
parameter plane (17, V), where 

^-(f^) = -^. (4.7) 

Localized solutions may only exist above the threshold V > Vthr{^) for any G M. The asymptotic 
approximation of the threshold (14. 7 1 is shown on Figure 4 by dotted curve. The threshold (I4.7t is an 
asymptotic approximation of the exact curve (14.31) . 

• When 7 = 0, the third-order ODE ( 14. 6t is related to the integrable Hirota equation IH73L which 
admits the exact single-humped localized solution = sech(Z). The single-humped solution 
exists everywhere on the two-parameter plane (17, V) above the threshold ( I4.7t . 

• When 7 = 1, the third-order ODE (14. 6t is related to the integrable Sasa-Satsuma equation IISS91I . 
which also admits the exact locahzed solutions everywhere on the two-parameter plane {0.,V) 
above the threshold (I4.7t . (The exact solution and its properties are given in Sections 4-5 of 
IISS91I .) The localized solution has a single-humped profile for Vthr(f^) < V < V|iumps(^^) and a 
double-humped profile for V > Vhumps(^) and 17 7^ 0, where 

(317)2/3 

Humps(J^) = (4-8) 



The threshold (14.81) in the ODE (14.61) corresponds to the condition A^ = A;^. It is shown on Figure 
4 by dashed-dotted curve. When 17 ^ and V > 0, the distance between two humps diverge and 
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the ODE has again the single-humped solution for = and V > with = ■^sech(Z). 
This solution corresponds to /c = in the ODE ( I4.6I) . 

• When 7 G M, the third-order ODE (I4.5t has the exact single-humped localized solution for = 
and V > with ^ = ■y=^sech(Z). The exact solution matches the member of the family of 
exact solutions in the two integrable cases 7 = and 7 = 1. It is shown in IPY05II by numerical 
analysis of the kernel of the linearization operator (see Sections 3-4 in IPYOSI ) that the family of 
single-humped localized solutions for 7 G K\{0, 1} is isolated from other families of localized 
solutions, i.e. the one-parameter family with 17 = and V > can not be continued in the 
two-parameter plane (0, V) as a single-humped localized solution. 




Figure 4: The boundary (I4.3t of the domain of existence for traveling solutions of the differential 
advance-delay equation (I1.12t and asymptotic approximations for the threshold on existence (14.71) (dot- 
ted curve) and the threshold on two-humped solutions (14. 8t (dashed-dotted curve). 

Results of the third-order ODE (14.51) give only necessary condition on existence of traveling single- 
humped solution in the differential advance-delay equation (I1.12t near the special values u = (vr — 2) //i^ 
and c = l/h, i.e. if the smooth solution (j){z) to the differential advance-delay equation (I1.12t exists, 
then it matches the analytical solutions of the third-order ODE i4.5l . Persistence proof of true localized 
solutions is a delicate problem of rigorous analysis, which is left open even for traveling kinks of the 
discrete 0^ equation IIP06I . (For comparison, persistence of solutions with exponentially small non- 
localized oscillatory tails was rigorously proved near the normal form equation in | IP06j .) 

When the discrete NLS equation (11.31) is the integrable AL lattice (11.71) . the exact traveling solutions 
exist in the form ^{Z) = sinh/t sech(KZ) for any k > and /5 € [0, 2tt]. These solutions match the 



20 



exact solutions for the Hirota equation for V > Vthv{^) and Q € M in the form ^'(C) = sech(C) with the 
correspondence k = eX. Additionally, when the discrete NLS equation (ll.3t has the nonlinear function 

/ = (1 - a)\Un\'^iUn+l + Un-l) + au^(u„+i + Un-l), a ^ 0, 

the exact traveling solution exists for k > and /3 = | in the form 

<I>(Z) = - sinh KPrh(K.Z). 

VI — 2a 

This solution matches the exact solution of the third-order ODE (I4.5t foiQ = and F > in the form 
"i!{C,) = —^^=sech{Z) (where 7 = —2as and a = 03). We do not know if there are any other examples 
of the discrete NLS equation il.3l with the nonlinear function ilAl which admit families of single- 
humped traveling solutions which match the analytical solutions of the third-order ODE (I4.5t . See 
IBGH04I for symbolic tests of explicit tanh and sech solutions of differential advance-delay equations. 



5 Conclusion 



We have considered stationary and traveling solutions of the discrete NLS equation (ll.3t with the ten- 
parameter cubic nonlinearity (II under the normalization constraint (ll.5t . We have proved in Lemma 
I2.3l that four constraints ( I2.7I ) are sufficient for reduction of the second-order difference equation (11.1 U 
to the first-order difference equation (I2.lt . The main result (Proposition B.lOt says that this reduction 
gives a sufficient condition for existence of translationally invariant stationary solutions that converge to 
the stationary solutions il.2l with c = in the continuum limit h ^ 0. 

Furthermore, we have proved in Lemma 13.11 that two additional constraints (13. 3t are sufficient for 
uniqueness of real-valued stationary locahzed solutions in the first-order difference equation (I2.lt . 
We have finally found in Section 4 that one more constraint gives the necessary condition for exis- 
tence of traveling solutions in the differential advance-delay equation ( I1.12t near the special values 
to = {tt — 2)/h? and c = \/h. Combining all constraints (I2.7t . (13. 3t . and (I4.4t and the normaliza- 
tion constraint ( ll.5t . we have a two-parameter family of the discrete NLS equation with the nonlinear 
function 

/ = (l-Q-2/3)|ii„|^(n„+i+ii„_i)+Qu^(n„+i+'u„_i)+/3(|n„+ip + |n„,_ip)(n„+i+n„_i), (5.1) 

where (a, /3) are two real-valued parameters. When a = /? = 0, the model (I5.lt is the AL lattice 
( ll.7t with the two-parameter family of traveling solutions in the domain (14. 3t . When a = 0, the model 
(I5.lt is related to the integrable Hirota equation (14. 5t with 7 = 0. When a = — i, the model (I5.lt is 
related to the integrable Sasa-Satsuma equation (I4.5t with 7 = 1. When /? = 0, the model (I5.lt admits 
one-parameter family of exact traveling solutions on the line u) = {irch — 2)//i^ and c > l/h. 



21 



Although the discrete NLS equation (ll.3t with the nonUnear function i5.ll has translationally invariant 
stationary solutions il.9i . the existence of traveling solutions (ll.lOt for any c 7^ is an open problem, 
which is left beyond the scopes of the present manuscript. Persistence of traveling solutions for small 
values of c can not be proved since infinitely many resonances with infinitely many Stokes constants 
appear in the limit c — > IIQPB05I . (It was shown in IQPBOSI with numerical computations of the 
leading Stokes constant that none of three particular discrete cj)'^ lattices exhibit families of traveling so- 
lutions bifurcating from the family of translationally invariant stationary solutions.) Although existence 
of translationally invariant stationary solutions gives only the necessary condition for persistence of trav- 
eling solutions, the integrable AL lattice represents at least one example of the discrete NLS equation 
(ll.7t where persistence of traveling solutions can in principle occur. 

Similarly, persistence of traveling solutions near the integrable normal form (I4.5t can not be proved 
since oscillatory tails are generic near the values uj = {tt — 2)/h? and c = l/h IIIP06II . While the AL 
lattice (15.11) with a = (3 = has exact traveling solutions between the two limiting cases c = and 
c = it is not clear if the translationally invariant NLS lattice (15.11) with a,/3 7^ exhibits any 
families of traveling solutions on the plane {uj,c). If such solutions exist for sufficiently small values 
of h, these solutions converge to the family (ll.2t in the continuum limit /i — > 0. This problem can be 
considered by means of numerical solutions of the differential advance-delay equation ( 11.121 ) similarly 
to the numerical works IAEHV05[fC06l . 
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